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Abstract 

Five physics mechanisms of interaction leading to binding of the molecular ion are identified, 
realized in a form of variational trial functions and their respective total energies are calculated. 
Each of them provides subsequently the most accurate approximation for the Born-Oppenheimer 
(BO) ground state energy among (two-three-seven)-parametric trial functions being correspond- 
ingly, H2-molecule plus proton (two variational parameters), H^-ion plus H-atom (three variational 
parameters) and generalized Guillemin-Zener (seven variational parameters). These trial functions 
are chosen following a criterion of physical adequacy, they include the electronic correlation in the 
exponential form ~ exp (77-12), where 7 is a variational parameter. Superpositions of two different 
mechanisms of binding are investigated and a particular one, which is a generalized Guillemin- 
Zener plus H2-molecule plus proton (ten variational parameters), provides the total energy at the 
equilibrium E = —1.3432 a.u. Superposition of three mechanisms: generalized Guillemin-Zener 
plus (H2 -molecule plus proton) plus (H^ -ion plus H) (fourteen parameters) leads to the total 
energy which deviates from the best known BO energy to ~ 0.0004 a.u. In general, our variational 
energy agrees in two-three-four significant digits with the most accurate results available at present 
as well as major expectation values. 
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I. INTRODUCTION 



We had learned recently from Takeshi Oka about a fundamental importance of the two- 
electron hydrogenic molecular ion H^~ in Physics, Chemistry and Astronomy [l| , in particular, 
that is a system of the lowest total energy among ones made from protons and electrons. 
The ion Hjj~ is one of the most abundant chemical compounds in the Universe being a 
major proton donor in chemical reactions in interstenar space. Experiment^, the H+ was 
discovered by J.J. Thomson in 1912 |2J. From another side, being theoretical physicists 
we are not aware about any popular quantum-mechanical textbook as well as a book on 
spectral theory where the molecular ion is even mentioned. Even though the H^" as the 
simplest three-center problem plays the same role of a fundamental object as the hydrogen 
and helium atom being the simplest one-center problems, and H^ and H2 being the simplest 
two-center problems. Needless to say that such a role implies to use Hjj~ as a test ground for 
different theoretical approaches. 

In general, the Coulomb system (3p, 2e) is very difficult for theoretical studies. Many 
theoretical methods were developed to study low-lying quantum states of this system. In 
particular, it became clear very quickly that interelectron correlation is of great importance 
and it has to be included into the variational trial function explicitly to assure a faster 
convergence. This conclusion is similar to one drawn by Hylleraas {^J for the He atom, and 
by James and Coolidge |4| for the H 2 molecule. Usually, the interelectron correlation was 
written in a form of monomial r™ 2 (Hylleraas jsj - James- Coolidge jj] form) or exp(— arf 2 ) 
(Gaussian form, see e.g. Ref. ja]). Recently, Korobov jg] showed in an explicit way that 
for the case of the Helium atom the use of the exponential form exp(— 7r 12 ) dramatically 
improves convergence and leads, in fact, to the most accurate results for the ground state 
energy for the Helium atom at present. Later on, it was shown that the similar use of 
the exponential correlation exp(— 7r 12 ) for the H 2 molecule allows us to construct the most 
accurate trial function among the few-parametric trial functions Q. A hint of why this 
ri 2 -dependence leads to the fast convergent results was given in 

In 2006, the Royal Society discussion meeting on H^ took place in London, UK (see [lj]) 
and at 2012 the Royal Society Theo Murphy Meeting on H^ held where different properties 
of the Hjj" ion and, in particular, various theoretical approaches to a study of the ion 
were presented (see |8[). It is worth mentioning that the benchmark results for the potential 



2 



energy surface of were reported by Adamowicz - Pavanello 9| . 

The goal of this contribution is twofold: (i) to study a physics of binding, and based on 
that (ii) to propose a maximally simple, compact and easy-to-handle trial function with few 
(non) linear parameters which would lead to a highly accurate Born-Oppenheimer ground 
state energy and major expectation values. We always assume that two-Coulomb charge 
effective interaction in a presence of other Coulomb charges is modeled by e~ ar , where a is 
a parameter (which is a lowest Coulomb orbital for the case of opposite charges), since a 
corresponding potential reproduces a Coulomb singularity at small distances and vanishes at 
large distances. Perhaps, it is worth noting that we are not aware of previous studies of the 
ion with trial functions involving r 12 in an exponential form with a single exception jic| 
where the Hjj" in a linear configuration in a magnetic field was explored. This paper can be 
considered as a continuation of a study of different physics mechanisms of binding started at 



where two possible mechanisms of binding and their superposition were explored (not 
dominant ones). 

In this paper atomic units (h — e — m e — 1) are used throughout. 

II. THE H+ ION IN THE BORN-OPPENHEIMER APPROXIMATION 

The Hamiltonian which describes the ion Hjj" under the assumption that the protons 
are infinitely massive (the Born-Oppenheimer approximation of zero order) and located at 
the vertices of an equilateral triangle of side R (see Fig. [TJ for the geometrical setting and 
notations), is written as follows: 

2 A2 



n 



2 r iA r 12 R 

J = l 3 = 1,2 J ' 

k=A,B,C 



where f)j = —iVj is the 3-vector of the momentum of the jth electron, the index k runs 
over protons A, B and C, r JjK is the distance between the jth electron and the nth proton, 
ri2 is the interelectron distance, and R is the interproton distance. 

It is a well established fact that the ground state of the Hjj" molecular ion is l 1 ^, an elec- 
tronic spin-singlet state, with the three protons forming an equilateral triangle in the totally 



symmetric representation A[ of a D 3 h point symmetry 12]. Thus, the ground state elec- 
tronic wavefunction should be symmetric under permutations of the three indistinguishable 
protons. This ground state is the major focus of the present study. 




FIG. 1: Geometrical setting for the hydrogen molecular ion H3 in equilateral triangular configu- 
ration. The three protons are located on the x-y plane forming an equilateral triangle with the 
origin of coordinates located at the geometrical center (circumcenter) of the triangle. 



It is worth mentioning that the best theoretical value at the moment for the Born- 
Oppenheimer ground state energy is E = —1.34383562502 a. u. [13] obtained with a basis 
of 1000 explicitly correlated spherical Gaussian functions with shifted centers. This value 
surpasses the previous record E = —1.343835624 a.u. by Cencek et al. which was obtained 
by using explicitly correlated Gaussian functions [li 



III. VARIATIONAL METHOD 



We use the variational procedure to explore the problem. The recipe for choosing the 
trial function is based on arguments of physical relevance, e.g. the trial function should 



support the symmetries of the system, has to reproduce the Cou 



omb singularities and the 



asymptotic behavior at large distances adequately (see, e.g. 15l-ll7l|). To the problem finding 



a bound state of (3p, 2e) giving rise to molecular ion we follow the description given in 



111- 

In practice, the use of such trial functions based on physical adequacy implies the conver- 
gence of a special form of the perturbation theory where the variational energy is the sum of 
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the first two terms of the perturbation series. Let us recall the essentials of this perturbation 



theory (for details, see 15Nl7l|). Let us assume that the original Hamiltonian has a form 
H = —A + V, where A is the Laplacian. As a first step we choose a trial function ?/;(*™ a ') 
(normalized to one) and find a potential for which such a trial function ^C" aZ ) is an exact 
eigenf unction, i.e. Vtriai — A^ tna ^ / ip^" 0,1 ' , with energy E tria i = 0. In a pure formal way 
we can construct a Hamiltonian Htriai = — A + Vtriai such that r H.tr%ai^ tTial ' ) =0. It can be 
easily shown that the variational energy 

E var = {i) itrial) \H\i) {trial) ) 

is nothing but the first two terms in the perturbation theory where the unperturbed problem 
is given by Htriai and the perturbation is the deviation of the original potential V from the 
trial potential V tr i a h namely, V pertU rbation = V — Vtriai- Eventually, we arrive at the formula 

E va r ^trial E\ {VperturbatioTi) > (2) 

here Ei(V pe rturbation) = (^ {tnal) \V P erturbation\'ip { - trtal) ) is the first energy correction in the per- 
turbation theory, where the unperturbed potential is Vtriai- It is worth noting that if the 
trial function is the Hartree-Fock function, the resulting perturbation theory is nothing but 
the Moeller-Plesset perturbation theory (see, e.g. jl9| . Section 15. 18) 1 . 

One of the criteria of convergence of the perturbation theory in V per t U rbation = V — Vtriai 
is a requirement that the ratio \ V per t U rbation/V\ should not grow when r tends to infinity in 
any direction. If this ratio is bounded by a constant it should be less than one. In fact, it is 
a condition that the perturbation potential is subordinate with respect to the unperturbed 



potential. The value of this constant contro 



s the rate of convergence - a smaller value of 



this constant leads to faster convergence 16j. Hence, the above condition underlines the 



importance of the large-range behavior of the trial functions. In the physics language the 
above requirement means that the phenomenon of the Dyson's instability should not occur 



(for a discussion see [15]) 2 . 



1 It is worth noting that the question about convergence of the Moeller-Plesset perturbation theory is not 



settled yet [20 1 

2 It is worth noting that this procedure for a selection of the trial function was applied successfully to a 
study of one-two-electron molecular systems in a magnetic field leading to highly accurate results. Many 



of these results are the most accurate at present (see 17] and 18[)- 
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In order to make the above-mentioned criteria of convergence concrete for the problem of 
Hjj~ let us introduce the exponential representation of the spin-singlet ground state function 

Xj/^) = (!-)- P 12 ) e -¥{r±A,r 1 B,r 1 c,r 2 A,r2B,r2C,yri2) ^ 

permjA, B,C} 

where ip(r 1 ,r 2 ) = v 9 ( r iA, r iB, r ic, r 2A, r 2B, r ic-, l r \2) is unknown phase. If d3J) is substituted 
to the Schroedinger equation a non-linear equation of a type of the multidimensional Riccati 
equation will occur. The analysis of the behavior of tp at large distances along any direction 
in six-dimensional physical space R 6 leads to the asymptotics of the form 

cp — > ar + O(logr) , r — > oo , (4) 

where the constant a can depend on the chosen direction, while the behavior near a Coulomb 
singularity r c of the potential is 

<p — > br + 0{r 2 ) , r — > rc , (5) 

where b makes sense of a cusp parameter (residue at the Coulomb singularity). Hence, 
a trial function with a phase Atrial which grows linearly at large distances guarantees a 
convergence of the perturbation theory where the chosen trial function is considered as zero 
approximation. Furthermore, for such a choice the ratio \(ftHai/ '^Pexact\ is bounded. Actually, 
this property implies that relative accuracies in energy and in an expectation value may not 
be very much different. 

There are the different ways to include electronic correlation explicitly in the trial wave 



function for two-electron problems. We mention three major approaches (see e.g. 21] ): the 
linear in ru, the gaussian exp(— arl 2 ) and exponential exp(7ri2) terms. Among them, the 
only factor exp(7ri 2 ) fulfills the physical adequacy requirements for a trial function described 
above and may lead to a convergent perturbation theory. 



IV. PHYSICS MECHANISMS OF BINDING 



As an illustration to what follows let us describe physics mechanisms of binding of the 
simplest two-center, one-electron molecular system H^. There are three of them: (i) a co- 
herent interaction of electron with two centers, (ii) incoherent interaction of electron with 
two centers, and (iii) a general, screened interaction of electron with two centers. A simplest 
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realization of these three mechanisms leads to the celebrated Heitler-London function (the 
product of two lowest, equally screened Coulomb orbitals, the electron "sees" two centers 
simultaneously), Hund-Mulliken function (the sum of two lowest, equally screened Coulomb 
orbitals, the electron "sees" one center only, it corresponds to the interaction H-atom with 
proton), and Gulliemin-Zener function (the symmetrized product of two lowest, unequally 
screened Coulomb orbitals, the electron "sees" both centers but differently). A linear super- 
position of these three functions gl ves the accuracy of order of 10- for any separation (see 
e.g. [121). 

Following the guidelines of Section III and the convergence requirement of the pertur- 
bation theory, the most general trial function for the ground state made out of individual 
exponentials is of the form: 



see Fig.l for notations, where the sum runs over the permutations of the identical protons 
A, B, C (S3 symmetry), and P\ 2 is the operator which interchanges electrons (1 -H- 2). 
The variational parameters consist of non-linear parameters «i_6 and 7 which characterize 
the ( ant i) screening of the Coulomb charges. The interproton distance R, see Fig.l is kept 
fixed and all there distances are chosen to be equal, thus, protons form equilateral triangle, 
R = R eq = 1.65 a.u. [10-11] 3 . The function (jSJ) is a symmetrized product of Is Slater type 
orbitals multiplied by the exponential correlation factor e 7ri2 . 

Calculations of the variational energy were performed using the minimization package M I- 
NUIT from CERN-LIB. Six-dimensional integrals which appear in the functional of energy 
were calculated numerically using a "state-of-the-art" dynamical partitioning procedure. 
The calculations were made in double-cylindrical coordinate system: (xi, yx, Zi,x 2 , 1/2, 22) — > 
(pi, 0i, z±, P2, 02, 22). As the first step the infinite domain of integration was reduced to a 
finite, compact domain in a form of two 3D cylinders in such a way that a contribution from 
the complement should be two orders of magnitude smaller than the requested absolute 
accuracy of integration. As the second step the domain of integration was subdivided into 

3 Almost all previous calculations were carried out for the same value of R eq . It is justified by the fact that 
the minimum of the potential curve in R is very flat being characterized by a small curvature. Hence, 
for a given accuracy in energy the accuracy of localizing the position of the minimum is much lower (see 





— air 1A — a2r 1B -a3r lc .— «4r'2A-"5»*2B-a6»~2C+7 r 'i2 



(6) 



perm{A,B,C} 



Table I) 



7 



972 up to 6534 subdomains following the profile of the integrand, in particular, separating 
out the domains with sharp changes of the integrand, the domains of the large gradients. 
Then each subdomain was integrated separately in a parallel manner with controlled abso- 
lute/relative accuracy (for details, see e.g. |17|). A realization of the routine required a lot of 
attention and care. During the minimization process both the domain and partitioning were 
permanently controlled and adjusted. Numerical integration of every subdomain was done 
with a relative accuracy of ~ 10 -3 — 10~ 7 depending on its complexity and relative contri- 
bution using an adaptive routine based on an algorithm by Genz and Malik [22J. This is an 
adaptive multidimensional integration routine (CUBATURE) of vector-valued integrands 
over hypercubes written by Steven G. Johnson with a vectorized prototype prepared by 



Dmitry Turbiner, http : // ab-initio .mit . edu / wiki / index. php / Cubature . Parallelization was 
implemented using the MPI library MPICH. Computations were performed on a Linux clus- 
ter with 48 Xeon processors at 2.67 GHz each, and 12Gb total RAM plus an extra processor 
serving as the master node. The complete minimization process for every Ansatz took from 
100 (single Ansatz, Eqs. (6-11), see below)) up to 2000 (triple mechanism Ansatz, Eq. (16), 
see below) hours of aggregated wall clock time. A single integration was about from 1 up to 
15 minutes. With optimal values of parameters, it took about 60 minutes (wall clock time) 
to compute a variational energy (calculating two integrals) with relative accuracy 10~ 7 . 



A. Single mechanism 

A goal of subsection is to identify some possible single mechanisms of interactions which 
occur in the system of three static protons and two electrons, (3p2e). 



1. Coherent interaction 



One of the most natural mechanisms of interaction is when each electron "sees" all charged 
centers simultaneously, coherently interacting with them. Such a mechanism has to be 
dominant at small internuclear distances in comparison with a characteristic length of the 
electron-proton interaction which seems to be of the order of 1 a.u. = 1 Bohr radius. A 
simplest realization of this mechanism in a wavefunction is given by a product of Coulomb 
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Trial Function 


E (Ry) 




a 2 


a 3 


Q'l 


a 5 


a 6 


7 


1. 


Coherent 


-2.583 


0.5711 


0.5711 


0.5711 0.5711 


0.5711 


0.5711 


0.29783 


2. 


Incoherent 


-2.595 


1.3557 











1.3557 





-0.06136 


3. 


H~ + 2p 


-2.148 


1.3284 








1.3284 








0.25244 


1 


H 2 +P 


-2.634 


0.9422 


0.9422 





0.9422 


0.9422 





0.43977 


5. 


H+ + H 


-2.650 


0.9114 


0.9114 











1.2499 


0.04359 


6. 


Generic 


-2.680 7 


-0.00353 0.18548 1.4245 1.0471 0.15082 0.58912 


0.21632 



TABLE I: The ground state energy of the (3p2e) system in different mechanisms of binding realized 
by Ansatz 1-6. 



orbitals 



coh 



-ai(r 1 A+r 1 B+ric+r2A+r2B+r2c)+'yri2 



(7) 



(Ansatz 1) which is a generalization of the Heitler-London function for the case of the 
molecular ion or the H2 molecule; it contains two variational parameters 011,7 both 
of them having a meaning of charge (anti) screening. Variational energy as well as the 
corresponding values of the parameters are presented in Table [B It is already a striking fact 
that even at R eq = 1.65 a.u. which is essentially larger than the Bohr radius the function (JJJ) 
leads to a binding of the (3p2e) system with sufficiently large ionization energy ~ 0.12 a.u. 
(H+^H 2 +p) 4 . 

2. Incoherent interaction 

Another natural mechanism of interaction is when each electron "sees" a single charge 
center, thus, realizing an incoherent interaction. In fact, it corresponds to the interaction 
of two hydrogen atoms and proton, H + H + p. Such a mechanism has to be dominant at 
large internuclear distances in comparison to a characteristic length of the electron-proton 
interaction of the order of 1 a.u. = 1 Bohr radius. A simplest realization of this mechanism 
is given by 

Ancoh = (I + P12) e- Ql(riA+r2fl)+ ^ ri2 , (8) 

peim{A, B,C} 



In order to calculate the ionization energy we use the H2 total energy -1.1745 a.u. (rounded) found in [2c 
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(Ansatz 2) which is a generalization of the Hund-Mulliken function for ion or H2 molecule, 
it contains two variational parameters 011,7 having a meaning of charge (anti) screening. 
Variational energy as well as the corresponding values of the non-linear parameters «i, 7 are 
presented in Table [B At R eq = 1.65 a.u. which is essentially larger than the Bohr radius, the 
function (jSJ) leads to a binding of the (3p2e) system with sufficiently large ionization energy 
~ 0.12 a.u. 4 . The parameter 7 takes a negative value indicating the effective attraction of 
the electrons in this mechanism. 

3. (almost) Incoherent interaction 

This mechanism of interaction appears when two electrons " see" the same charge center 
interacting coherently, then, the interaction with two other charged centers appears as a 
result of the exchange interaction. In fact, it corresponds to the interaction of negative 
hydrogen ion with two protons, H~ + 2p. Such a mechanism has to be important at large 
internuclear distances comparing to a characteristic length of the electron-proton interaction 
of the order of 1 a.u. = 1 Bohr radius. It has to be suppressed due to a small probability 
of the coherent interaction of two electrons with same proton. A simplest realization of this 
mechanism is given by 

ip^ncoh = (I + P12) e- a *"^+™ , (9) 

perm{A,B,C} 

(we call it the Ansatz 3). Perhaps, it is not surprising that such a trial function does not 
lead to a binding (see Table I). 

4- H2-molecule + proton interaction 

In this case both electrons "see" two protons forming H 2 molecule which interacts with 
proton. Though the interproton distance R eq = 1.65 a.u. for Hg is larger than the equilib- 
rium distance for the H2 molecule, which is equal to ~ 1.4 a.u., it can be described by the 
Heitler-London function for the H2 molecule. Eventually, as a realization of this mechanism 
we choose 

^ H2 = (1 + p 12 ) e - ai{riA+riB+r2A+r2B)+ ^ ri2 , (10) 

perm{,4,B,C} 
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(Ansatz 4) which contains two variational parameters «i,7 having a meaning of charge 
(anti)screening. Variational energy as well as the corresponding values of the non-linear pa- 
rameters ai, 7 are presented in Table [B The function (1101) leads to a binding with sufficiently 
large ionization energy ~ 0.15 a.u. It is worth noting that an attempt to improve a descrip- 
tion of the H 2 molecule by replacing the Heitler-London function by the Guillemin-Zener 
one 

e -a 1 (r 1A +riB+r2A+r2B) e _a M ( r iA+riB)-ai,2( r 2A+r 2 B) 

leads to unessential improvement in ionization energy ~ 0.16 a.u. 4 . 

5. H^-ion + H-atom interaction 

In this case one electron "sees" two protons forming H^ ion and the second electron 
"sees" the third proton forming H-atom. The interproton distance R eq = 1.65 a.u. for 
is smaller than the equilibrium distance for the H^ ion, which is equal to ~ 2. a.u., H^ 
contribution can be described by the Heitler-London function for the H^ ion. Eventually, 
as a realization of this mechanism we choose 

^ H + H = (1 + P 12 ) e - ai(riA+ris) - a6r2C+7ri2 , (11) 

pevm{A,B,C} 

(Ansatz 5) which contains three variational parameters a\, a^, 7 having a meaning of charge 
(anti)screening. This function can be considered as a generalization of the celebrated 
Guillemin-Zener function written for the H^-ion and the H 2 -molecule. Variational energy 
as well as the corresponding values of the parameters are presented in Table HI The func- 
tion (TTTl) leads to a binding with sufficiently large ionization energy ~ 0.14 a.u. 4 which is 
nevertheless smaller than one for the Ansatz 4- The parameter 7 is close to zero indicating 
effectively the small repulsion of the electrons. 

6. Generic interaction 

Naturally, a general mechanism of interaction is related to interaction of an electron with 
a proton (or another electron) with its particular (anti) screening due to the presence of other 
charges. A simplest realization of this mechanism is given by the function ([6]) (see above). 
The variational parameters consist of non-linear parameters ai_6 and 7 which characterize 
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the (anti)screening of the Coulomb charges. Appropriate degeneration of these parameters 
allow us to reproduce different mechanisms of interaction 1.-5. Variational energy as well as 
the corresponding values of the parameters are presented in Table HI The function ([H]) leads 
to a binding with quite high ionization energy ~ 0.17 a.u. 4 . 



B. Mixed mechanisms 



In previous subsection we identified possible single mechanisms of interactions which oc- 
cur in the system of three static protons and two electrons, (3p2e). Based on the values 
of ionization (binding) energies these mechanisms can be classified coming from the highest 
ionization energy at Generic interaction mechanism (Ansatz 6.) followed by the mecha- 
nism Hz -ion + H-atom interaction (Ansatz 5.) and then Hz-molecule + proton interaction 
(Ansatz 4-) ■ The goal of this subsection is to consider different superpositions of single 
mechanisms. The mechanism (almost) Incoherent interaction (H~ + 2p, Ansatz 3.) is ex- 
cluded as one giving no binding to the (3p2e) system. In all mixed mechanisms the Generic 
interaction mechanism will be always present as one giving the highest ionization energy 
among single mechanisms. 



1. Generic interaction plus coherent interaction 



Let us consider a superposition of the Generic interaction (Ansatz 6.) and Coherent 
interaction (Ansatz 1.) (see 



Oil 



4>g + A coh ip coh . 



(12) 



Formally, the trial function depends on ten parameters (one linear and 9 non-linear). 
The results are presented in Table [TTJ It is interesting to note that the non- linear parameters 
are almost unchanged in comparison to the parameters of the individual Ansatze 1. and 6. 
with a single exception of the parameter 7. If it is compared with the Ansatz 6. this Ansatz 
increases the ionization energy to ~ 0.001 a.u. 
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E (Ry) 


a.\ Q2 «3 «4 Q-5 «6 7 


A a 7 


Generic Ansatz (6) 


-2.680 7 


-0.00353 0.18548 1.4245 1.0471 0.15082 0.58912 0.21632 




(6) e (7) 


-2.683 2 


-0.00294 0.21022 1.3849 1.0199 0.17103 0.59084 0.26044 


-0.51154 0.59589 0.86229 


(6) © (8) 


-2.683 5 


-0.02840 0.19525 1.4780 1.06860 0.12045 0.52108 0.24008 


-0.03987 1.71230 0.69991 


(6) © (11) (H+ + H) 


-2.685 4 


-0.05699 0.22600 1.40290 1.04611 0.17386 0.55836 0.14867 


-1.52370 1.50704 -1.21701 
1.10855 


(6) © (10) (H 2 ) 


-2.686 4 


0.00515 0.19795 1.4043 1.03200 0.13466 0.59350 0.13686 


-2.34310 0.99547 -0.46046 


(6) © (10) (H 2 ) © (11) (H+ + H) 


-2.686 852 


0.00040 0.23034 1.42450 1.00418 0.12283 0.56933 0.13409 


-2.34400 0.94900 -0.55946 
0.16243 1.72960 -0.87202 
0.066457 


Best value [13] (rounded) 


-2.687 671 







TABLE II: The ground state energy of HJ at R eq = 1.65 a.u. and the non-linear variational parameters in [a.u.]^ 1 corresponding to the trial 
function ([6]) and different degenerations of it. For the case (6) + + H, the a parameter in the second row corresponds to the effective 
charge for the H-atom part. 



2. Generic interaction plus incoherent interaction 

Let us consider a superposition of the Generic interaction (Ansatz 6.) and Incoherent 
interaction (Ansatz 2.) 

1p62 = IpG + Ancohlpincoh ■ (13) 

Formally, the trial function i/jq 2 depends on ten parameters (one linear and 9 non-linear). 
The results are presented in Table [Til It is interesting to note that the non- linear parameters 
of the Ansatz 6. in ffl3|) are almost unchanged in comparison to the parameters of the 
individual Ansatz 6. with a single exception of the parameter a% while the parameters of 
the Ansatz 2. are changed significantly. If it is compared with the Ansatz 6. this Ansatz 
increases the ionization energy to ~ 0.002 a.u. 

3. Generic interaction plus H2-molecule + proton interaction 

Let us consider a superposition of the Generic interaction (Ansatz 6.) and H 2 -molecule 
+ proton interaction (Ansatz 4-) 

^64 = 4>g + A H2 ?p H2 . (14) 

Formally, the trial function t/> 64 depends on ten parameters (one linear and 9 non-linear). 
The results are presented in Table (TTJ It is interesting to note that the non- linear parameters 
of the Ansatz 6. in (fT4"j) are almost unchanged in comparison to the parameters of the 
individual Ansatz 6. with a single exception of the parameter 7 while the parameters of 
the Ansatz 3. are changed insignificantly. If it is compared with the Ansatz 6. this Ansatz 
increases the ionization energy to ~ 0.003 a.u. In comparison with the most accurate total 
energy this Ansatz reproduces 3 figures. 

4. Generic interaction plus H^-ion + H-atom interaction 

Let us consider a superposition of the Generic interaction (Ansatz 6.) and H^-ion + 
H-atom interaction (Ansatz 5.) 

^65 = IpG + A H + +H lP H + +H . (15) 
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Formally, the trial function ^5 depends on eleven parameters (one linear and 10 non-linear). 
The results are presented in Table [Til It is interesting to note that the non-linear parameters 
of the Ansatz 6. in (1151) are almost unchanged in comparison to the parameters of the 
individual Ansatz 6. with a single exception of the parameter 7 while the parameters of the 
Ansatz 4- ar e changed significantly, especially the parameter 7. If it is compared with the 
Ansatz 6. this Ansatz increases the ionization energy to ~ 0.002 a.u. In comparison with 
the most accurate total energy this Ansatz reproduces 3 figures. 

5. Generic interaction plus H2-molecule + proton interaction plus H^-ion + H-atom interac- 
tion 

Following the analysis of the variational total energies found in single mechanisms of 
binding of the H^"-ion we arrive at the hierarchy: Generic interaction (Ansatz 6.) goes 
first, then H^-ion + H-atom interaction (Ansatz 5.) is the second and H.2-molecule + 
proton interaction (Ansatz 4-) as the third. In turn, the analysis of the variational total 
energies found in mixed mechanisms of binding involving two mechanisms we arrive again 
at a conclusion that above-mentioned three mechanisms of binding play a dominant role. 
Thus, a natural move is to consider a superposition of these three mechanisms of binding: the 
Generic interaction (Ansatz 6.), H^-ion + H-atom interaction (Ansatz 5.) and H 2 -molecule 
+ proton interaction (Ansatz 4-), 

^645 = + A H2 tjj H2 + A H + +H 1p H + +H • (16) 

The trial function ip^^ depends on fourteen parameters (two linear and 12 non-linear). The 
results are presented in Table (TTJ It is interesting to note that remarkably the non-linear 
parameters of the Ansatz 6. in (fTBT) are almost unchanged in comparison to the parameters 
of the individual Ansatz 6. itself as well as in Ansatz (fi~5l) and Ansatz (fTij) with a single 
exception of the parameter 7 which nevertheless remains almost unchanged for Ansatz (fi~5l) , 
Ansatz ( fl^l) and Ansatz ( TTBl) . The non-linear parameters of the Ansatz 5. are slightly 
changed in the Ansatz ( [151) and Ansatz ( Tl6l) . We consider this as an indication that the 
Generic interaction is described adequately by the Ansatz 6. ([6]) as well as the H 2 -molecule + 
proton interaction is described adequately by the Ansatz 4- ([10]) . The non-linear parameters 
of the Ansatz 5. are changed significantly, especially the parameter 7, comparing to ones in 
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the Ansatz (TH|) and Ansatz (fl~6]) . Also the non-linear parameters of the Ansatz 4- in the 
Ansatz (j!4p and Ansatz ffl6|) are very much different. We consider this as an indication that 
the mechanism of the H^-ion + H interaction is described inadequately by the Ansatz 4- 
( TTTj) . Perhaps, for H^-ion part the Heitler- London function for H^" should be replaced by 
more adequate Guillemin-Zener function 

/ \ (1) (2) 

e -ct4,i{riA+riB) e - a \A r iA-a\ [r 1B 

(see (ITT]) ). It will be done elsewhere. 

If the Ansatz f[T6"j) is compared with the Ansatz 6. this Ansatz leads to increase in the 
ionization energy to ~ 0.003 a.u. Making a comparison of the obtained variational total 
energy with the most accurate total energy it appears that this Ansatz provides a deviation 
from the best known energy which is equal to ~ 0.0004 a.u. (see Table IITT1) . 



V. RESULTS 

In Table ILTTI we present the results for the ground state energy at interproton equilibrium 
distance of the Hjj" molecular ion obtained previously using different methods. It can be seen 
that the Born-Oppenheimer ground state energy obtained using the trial function ( TT6|) is the 
most accurate (the lowest) energy obtained with a few parametric functions. In particular, 
the trial function (f!6l) gives a lower energy than the energies obtained with the explicitly 



correlated functions based on both Gaussians in r\% [24| and linear in r%2 |25|], when a 
relatively small number of terms with non-linear parameters is involved. The trial function 
(Fl6|) is more accurate than all(!) traditional CI calculations which were performed before 



1990 (see 26]) even including the one with the largest set of 100 configurations 23, 128] 5 
and as well as [32], however, being worse than the largest CI calculation ever performed 
[33] (700 configurations). In those CI calculations no explicit correlation was included. 
The trial function (fl~6]) is more accurate than some CI calculations with explicitly involved 
r 12 dependence [34] (36 configurations) and [35] (192 configurations) but gets worse when 
(much) larger number of configurations is included [36-37]. In general, the 14-parametric trial 
function ( fl~6l) gives the eighth best result among known in literature so far (see Table IIIII) 



5 for a list of 42 calculations of the ground state energy of in the period 1938-1992 see Rcf. 26], for a 
list of selected ab-initio calculations till 1995, see [5j 
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being incomparably simpler with respect to all known trial functions. It is worth noting 



that the 7-parametric trial function fl6]) 
in literature (see Table HTT1) . 



11] gives the thirteenth best result among known 



The list of major expectation values obtained using different trial functions ([6]), (12) - 
(16) and its comparison with results of other calculations is given in Table IIVI A reason- 
able agreement for the expectation values is observed. In particular, for the expectation 
values (1/tia), (x 2 ), (z 2 ) and (r 2 ) agreement within ~ 1% with ours and all other calcu- 
lations is observed, including ones obtained in the large CISD-R12 calculations 25) . Also, 
for the expectation value of (l/r 12 ) we have an agreement in the first significant digit with 
other calculations, and we are also in close agreement to the value obtained with the corre- 
lated Gaussian (unrestricted) wavefunction with 15 terms 24j. For the expectation values 
(1/tia) and (z 2 ) we observe agreement with other calculations in 3 and 2 significant digits, 
respectively. These facts seem to indicate that the presented expectation values are very 
accurate, corroborating the quality of the trial functions ([6]), (12)-(16) which gives 2-3 s.d. 
correctly. It is worth noting that since there are no criteria about accuracy of the obtained 
expectation values, we can only comment about the agreement within our results based on 
different trial functions and those obtained by other approaches. Since the exact phase of 
the ground state eigenfunction is uniformly approximated by the phases the trial functions 
(6), (12) - (16) we guess that the relative accuracy reached in variational energy has to be 
similar to the relative accuracy in an expectation value. It implies that the expectation 
value based on (6), (12) - (15) should provide two significant digits correctly, while (16) 
should give three significant digits (see Table lIVT) . Following this assumption we can state 
that < r 12 >= 1.9(9) a.u., < l/r 12 >= 0.62(6) (a.u.)^ 1 (hence, in Ref . [27] a single figure 
was found correctly), < 1/t\a >= 0.85(5) (a.u.)^ 1 (it agrees in 2 figures with all previous 
calculations, see Table HVl) . < x\ >= 0.76(1) (a.u.) 2 (it agrees in 2 figures with all previous 
calculations, see Table HVl) . < z\ >= 0.54(2) (a.u.) 2 (it agrees in 2 figures with all previous 
calculations except (12) and (14), see Table HVj) . < r\ >= 2.0(6) (a.u.) 2 (it agrees in 2 figures 
with all previous calculations, see Table HVT) 
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E (a.u.) 


R (a.u.) 


method 


reference 


-1.339 7 


1.66 


CI-GTO, > 120 configs 


[29] (1970) 


-1.306 29 


1.65 


GG, 3 terms, 5 non-linear params 


[241 (1973) 


-1.327 25 


1.65 


GG, 6 terms, 7 non-linear params 




-1.331 47 


1.65 


GG, 10 terms, 9 non-linear params 




-1.332 29 


1.65 


GG, 15 terms, 11 non-linear params 




-1.334 382 


1.65 


R12, 10s basis set 


[25] (1993) 


-1.334 632 


1.65 


R12, 30s basis set 




-1.340 34 


1.65 


7-Parametric Trial Function ([6]) 


[11] (2011) 


-1.340 5 


1.6405 


CI -GTO, 48 configs 


[27] (1971) 


-1.340 5 


1.65 


CI -STO, 100 configs 


[28] (1971) 


-1.342 72 


1.65041 


CI-GTO, 108 terms 


[31] (1985) 


-1.342 784 


1.6504 


CI-GTO, 8s3pld/[6s3pld] basis set 


[32] (1978) 


-1.343 40 


1.6504 


CI-GTO, Ws4p2d basis set 


[33] (1990) 


-1.343 822 


1.65 


CI-GTO, 700 terms 


[34] (1990) 


-1.342 03 


1.6504 


CI with rl2, 36 configs 


[35] (1982) 


-1.343 422 


1.6504 


CI with rl2, 192 configs 


[361 (1984) 


-1.343 426 


1.65 


14-parametric Trial function (1161) 


present 


-1.343 500 


1.6504 


CI with rl2, 13s3p/[10s2p] basis set 


[37] (1988) 


-1.343 828 


1.65 


CI with rl2, 13s5p3d basis set. 


[38J (1990) 


-1.343 835 


1.65 


R12, 30s20pl2ci9/ basis set 


[25] (1993) 


-1.343 35 


1.65 


GG, 15 terms, 135 non-linear params 


[24] (1973) 


-1.343 835 624 


1.65 


GG, 600 terms 


[14] (1995) 


-1.343 835 625 02 


1.65 


ECSG, 1000 terms 


[13] (2009) 



TABLE III: A selection of the calculations for the Born-Oppenheimer ground state energy at 
equilibrium distance of Hjj". Record calculations of Ref. [r| (2009) and Ref. \l\ (1995). CI denotes 
Configuration Interaction, STO - Slater Type Orbitals, GTO - Gaussian Type Orbitals, GG - 
correlated Gaussians (Gaussian Geminals), R12 - the CI calculation augmented by terms linear 
in ri2, ECSG - Explicitly Correlated Spherical Gaussian functions. Our present results: Trial 
Function (|16p with fourteen parameters optimized. 



18 



Expectation Value 


Trial Function 


Others 


/ \ 

\ r l2) 


2.0032 ((6J) 






1.9931 (H2D 






1.9927 CGD 






1.9831 (fT5l) 






l.yooo l|l4|) 






l.yoo4 (|lb|) 




( l l r l2) 


U.oolo (|6J) 


0.5941O 




0.6302 (HI 


0.62636 c 




0.6293 ((ISI) 






6273 (fT5l) 






n coon /i 1 /i i\ 
U.ozoU {(141) 






U.OZOl (|1D|) 






0.8548 Q6J) 


0.85519 




0.8549 (H2D 


0.8553 e 




0.8553 (fH 






0.8555 ([E]) 






0.8556 CED 






0.8553 (USD 





TABLE IV: Expectation values (in a.u.) for the Hjj~ ion in its ground state obtained with the trial 
functions © and (|12M16p . Corresponding results obtained with other methods are displayed for 
comparison. The coordinates x, y, z and r are measured from the center of the equilateral triangle 



formed by protons. " Ref.j^l CI- .13: '' CI wavefuncion (I) in He I'. 



Correlated Gaussian 



wavefunction with the Ws8p6d4:f basis set in Ref. 
VI. CONCLUSION 



25]. 



(unrestricted) wavefunction with 15 terms in Ref. 24 1; d CI wavefunction in Ref. 



30j; e CISD-R12 



We identified five mechanisms of interaction leading to a binding in the Coulomb problem 
of two electrons in the field of three fixed charged centers (2e3p) in equilibrium geometry 
- the charged centers form the equilateral triangle with side 1.65 a.u. Among these five 
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Expectation Value 


Trial Function 


Others 


/ 2\ 1 2\ 


0.7711 (|6J) 


0.75818 




0.7666 (H2D 


0.7595 b 




0.7610 CGD 


0.75913 c 




7588 (fl5l) 


75968 d 




n Vfino /i i 1 r\ 
U./DUy (|14|) 


U. < OUD 




U.rblo (|lb|) 




/-v2\ 
V> ) 


n conn /f7TI\ 

u.ooyy (|oJ) 


0.54o02 




0.5337 (HI 


0.5451 6 




0.5420 (d3I) 


0.54085° 




5375 (fl"5l) 


0.54179 d 




(J.5o4o {(141) 


0.5396 




U.04ZU (|lD|) 




/„2\ 
(r ) 


2.0822 Q6J) 


2.06442 




2.0669 (H2D 


2.0640 b 




2.0640 (HI 


2.05911 c 




2.0550 CED 


2.06114 d 




2.0565 (HID 






2.0646 (USD 





TABLE IV: Continued. Expectation values (in a.u.) for the Hjj~ ion in its ground state obtained 
with the trial functions ([6]) and (|12M16p . Corresponding results obtained with other methods 
are displayed for comparison. The coordinates x^z and r are measured from the center of 



the equilateral triangle formed by protons. a Ref. 



Ref. 



30]; e CISD-R12 wavefunction with the 10s8p6d4/ basis set in Ref. 



25]. 



2j CI-43; b CI wavefuncion (I) in Ref. 
Correlated Gaussian (unrestricted) wavefunction with 15 terms in Ref. 24]; d CI wavefunction in 



mechanisms three lead to a strongest binding: H 2 -molecule plus proton, H\ion plus H- 
atom and generic. Each mechanism of interaction is easily modeled by a certain two-three- 
seven-parametric wavefunction, respectively, which is taken as a variational trial function. 
Corresponding variational energies are certainly the lowest ones among two-three and seven 
parametric trial functions, respectively (see Table [TTT]) . 
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In particular, a simple and compact 7-parametric variational trial function fl6]) (generic 
mechanism) already provides surprisingly accurate Born-Oppenheimer energy for the ground 
state of such a complicated molecular system as fhj~. The minimum energy is found to be 
E = —1.34034 a. u. at an equilibrium interproton distance R = 1.65 a. u. This result for 
the energy is the most accurate among the values obtained with several parametric trial 



functions of GG type 24J (Gaussian in r^), when a relatively small number of terms and 
non-linear parameters are involved, see Table IIHI In particular, it is more accurate than the 



energies obtained with the explicitly correlated approaches of Ref. [25]. 



It seems natural to advance by considering a superposition of three dominant mechanisms 
of binding like it was done in [7] for H 2 -molecule, 

* = ip G + A H2 tp H2 + A H + +H ip H + +H , 

as well as a superposition of two dominant mechanisms of binding. It immediately gives an 
essential improvement in the energy (see Table HTl) . In particular, the function (16) (a super- 
position of three dominant mechanisms), which contains fourteen variational parameters, 



allows us 


to 


get 


Woznicki 


36 


in 


than one by 


37 1 



in CI with rn included using 192 configurations. Our result is slightly worse 



371 ] in CI with r 12 included, using 13s3p/[10s2p] basis set. 
The obtained energy (16) is among eight the most accurate variational results ever calcu- 
lated so far but using incomparably much more simple trial function. The major expectation 
values in Table [IV] are gradually changed with the moves from one Ansatz to another seem- 
ingly demonstrating a convergence. It is clearly seen that in any expectation value (as well 
as in BO energy) two significant digits remain stable, they do not depend on the Ansatz 
considered. Thus, it can be stated that two significant digits in any expectation value are 
calculated reliably. Without doubt, the trial function ([6]) can be used to study potential 
energy surface. It seems that two significant digits are guaranteed. As an illustration on 
Fig. 2 the potential curve E vs R is presented. It is based on a linear superposition of Ansatz 
© and ©. 

It is worth emphasizing that the main virtue of functions ([6]), ( TTBl) is their compact- 
ness. A similar analysis of binding can be carried out for low-lying states and corresponding 
wavefunctions written. In particular, the function ([6]) can be easily modified for a study 
of spin-triplet states, as well as the low-lying states with non-vanishing magnetic quantum 



21 



-1 



h- 
LU 



-1.05 



-1.1 



-1.15 



-1.2 



-1.25 



-1.3 



-1.35 



% (fixed) + 1- coh 
*F 1 4 params 
Meyer et al. 



i- 

LU 




1.25 1.5 1.75 2 2.25 2.5 2.75 
R (a.u.) 



4 

R (a.u.; 



FIG. 2: Potential Energy Curve for the ion in triangular equilateral configuration as a function 
of the triangle size R obtained with the 7+3 parameter Ansatz f 1 1 2 [) (see [111]). Results by Meyer 



et al. obtained using CI 10s, 4p, 2d GTO basis 
at equilibrium is marked by bullet 



331 ] are shown for comparison. The energy of (16) 



number. A generalization to more-than-two electron molecular systems seems also straight- 
forward. 

In order to conclude we must mention that a domain of applicability of BO approximation 
of the zero order (infinitely massive particles of positive charge are assumed) is limited. 
This question can be thorough studied on the example of two-center problem exploring 
the H^~, D^~, ions. We used the ground state energies and the expectation values for 
the internuclear distance obtained in the Lagrange-mesh method for two-center problem 
EBo,Req and the H;]", D^~, T^~ ions [39[ with relative accuracy ~ 10~ 12 (for energies) and 
~ 10~ 7 (for the internuclear distances) found the fits 



E(M, m e 



E 



BO 



0.468436 



ira 



-j^ , < R NN (M,m e ) > -R eq 



2.84964 J— , 



where m e is the electron mass and M is the mass of a heavy particle. Hence, the finite-mass 
corrections begin to contribute to the 2-3rd significant digit. One way to include the finite- 
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mass corrections is to introduce the zero energy state corresponding to the lowest energy 
state of nuclear motion. In this case the zero energy is equal to the BO ground state energy 
of the zero order (a minimum at the potential curve) plus the first vibrational energy. One 
of the ways to find the first vibrational state is to consider a harmonic approximation of the 
bottom of the potential curve. One can see that the zero energy coincides to the non-BO 
energy up to four-five significant digits. In order to get further improvement it is necessary 
to go beyond harmonic approximation, also taking into account non-adiabatic corrections 
and write more accurate BO trial functions. It makes the overall analysis complicated. We 
think it is more relevant to consider an alternative way treating the problem solving the 
exact 3-body Schrodinger equation without any type of BO approximation similar to what 



was done for example in 



41] . In particular, it would allow to study the essential effects 



which do not exist in BO approximations like the electric quadrupole transitions. 

Similar consideration can be performed for three-center problem finding Ebo, Reg an d for 



40] it can be seen that finite-mass 



the H3 , D3 , T3 ions. From non-BO calculations made in 
corrections are essentially larger than for ones for two-center case. They contribute to the 
2-3rd significant digits to the ground state energies of hLt , D3 , T3 and the expectation value 
for internuclear distance, 

E(M,m e ) - E BO = 1.27232 J , < R NN {M,m e ) > -R eq = 4.21669 , 

where m e is the electron mass and M is the mass of a heavy particle. Taking into account 
the lowest vibrational energy found in Q the agreement between zero energy and non-BO 



ground state energy 



40J is improved up to 4th significant digit only. We would consider it 



as indication that the exact, 5-body Schroedinger equation has to be treated when going 
beyond the fifth significant digit in energy. The consideration presented in |40j] shows that 
is feasible. It seems evident that the dominant mechanisms of interaction leading to binding 
found in BO approximation will remain dominant beyond of it in non-BO consideration. 
We think that namely in such an approach the radiative transitions in the ion have to 
be studied. 
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